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ABSTRACT 

The Mbssbauer effect has been used to investigate electronic 

shielding by closed electron shells in salts of trivalent thulium, by 

measuring the temperature dependence of the nuclear quadrupole 

1 9 

splitting of the 8.4Z keV gamma transition in Tm . The nuclear 
quadrupole interaction was studied for Tm^ ions in thulium ethyl 
sulfate, thulium oxide and thulium trifluoride within a temperature 
range from 9.6°K to 1970°K. The interpretation of the experimental 
data in terms of the cont r ibutions of distorted closed electron shells 
to the quadrupole interaction yields values for electronic shielding 
factors. The results lead to amounts of 10% or less for the atomic 
Stcrnheimer factor Rq. The experiments also reveal substantial 
shielding of the 4f electrons from the crystal electric field, ex - 
pressed by the shielding factor tr 7 . Values of Z50 and 1Z8 are ob- 
tained for the ratio ( 1 -y^) /( 1 -cr^) for thulium ethyl sulfate and thu- 
lium oxide respectively, where is the lattice Sternheimer factor. 
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I. INTRODUCTION 

The technique of recoilless nuclear resonance absorption of 
gamma radiation, the so-called MHssbauer effect (1), has been em- 
ployed in numerous experiments in recent years (Z), (3). By binding 
a radioactive nucleus in a crystal lattice the emitted gamma radiation 
will, under certain conditions (Z), have essentially the natural line 
width as determined by the Heisenberg uncertainty relation and an 
energy exactly equal to the excitation energy of the nucleus. If a 
nucleus of the same isotope which is in its ground state is also bound 
in a lattice, there is a large probability for nuclear resonance ab- 
sorption of the gamma radiation. A distinct advantage of this tech- 
nique is the inherently high energy resolution that is available. For 
example in the experiments to be described here, the resolution is 
one part in 10^. Energy resolutions of this order make it possible 
to study nuclear properties as well as solid state effects in the crys- 
tals that are used to bind the nuclei. We make use of the MBssbauer 
effect here to study the nuclear hyperfine interactions in salts of 
rare earths, specifically thulium salts. 

Measurements of the nuclear quadrupole interaction in salts of the 
rare earth elements yield information on the quadrupole moments of 
the relevant nuclear states and on the electric field gradients which 



(1) R. E. Mdssbauer, Z. Physik 151, 1 2 4 ( 1958); Naturwissen- 
schaften 45, 538 (1958); Z. NcTmTfor sch. 14a, Z 1 1 ( 1959) 



(Z) See for instance H. F rauenfe lde r , The MBssbauer Effect , 

(W. A. Benjamin Inc. , New York, 1 9&Z) 

(3) The Proceedings of the Third International Conference on the 

Mdssbauer Effect appear in Rev. Mod. Phys . 3_6, 333-504 ( 1964) 
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exist in the salts at the nuclear sites. The extraction of the compo- 
nents of the electric field gradient tensor from such measurements 
is rather straightforward if the values of the nuclear quadrupole mo- 
ments have been obtained by other methods such as Coulomb excita- 
tion techniques. On the other hand the determination of nuclear mo- 
ments of rare earth nuclei by measurements of the nuclear quadru- 
pole interaction is rather involved since this requires a calculation 
of the components of the electric field gradient tensor at the nuclear 
sites. A calculation of the electric field gradients for salts of the 
rare earths can be performed at present only with limited accuracy. 
Uncertainties in excess of 30% are typical. It therefore appears 
that measurements of the nuclear quadrupole interaction in solids of 
the rare earths are at present of more importance for studies of the 
sources of the electric field gradients than for determination of 
nuclear quadrupole moments. 

The electric field gradient at the nuclear site of a certain ion 
originates from a number of different sources. Major sources are 
distortions of the electronic shells of the ion. These distortions 
result from the interactions of the electrons of the ion with the crystal 
electric field (CEF) produced by the surrounding ions in the lattice, 
provided the arrangement of the surrounding ions reflects a point 
symmetry lower than cubic. The field gradient at the nuclear site 
results not only from the distorted partially filled 4f electron shell 
of the rare earth ion, but also from distorted closed electron shells. 
These distortions of the closed electron shells of the ion constitute a 
major source of uncertainty in calculations of the electric field 
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gradient at the nuclear site. The deviations of the closed shells from 
spherical symmetry (electric multipole polarization) usually lead to 
substantial reduction or enhancement (shielding or antishielding) of 
the electric field gradient at the nuclear site. Sternheimer (4), (5) 
was first to emphasize the importance of magnetic dipole and electric 
quadrupole polarizations of closed shells and pioneered in calculating 
the contributions of closed shell deformations to the nuclear hyper- 
fine interactions. 

The nuclear quadrupole interaction depends strongly on the elec- 
tronic state of the ion. The electronic states which arise when a 
rare earth ion is incorporated in a crystal lattice are basically 
caused by the interaction of the CEF and the electrons in the partially 
filled 4f electron shell, but the splittings of these electronic levels 
are also strongly influenced by distortions of the closed electron 
shells (6), (7), (8). In order to account for the modification of the 
CEF splitting which results from electronic shielding, one has to 
consider the quadrupole moment as well as higher multipole mo- 
ments induced in the closed shells. 



(4) R.M. Sternheimer, Phys . Rev. 80, 102(1950); 105 , 158 

( 1957); R.M. Sternheimer and H. M. Foley, ibid. 102 , 731 
(1956); H. M. Foley, R.M. Sternheimer and D. Tyko, ibid. 

92> *734 ( 1954) 

(5) R.M. Sternheimer, Phys. Rev. 84_, 244 ( 1951); 9 5 , 736 ( 1954) 

(6) D.T. Edmonds, Phys. Rev. Letters J_0, 129 (1963) 

(7) R. G. Barnes, E. Kankeloit, R. L. Mbssbauer and J. M. 
Poindexter, Phys. Rev. Letters 1 1 , 253 ( 1963) 

(8) J. Blok and D. A. Shirley (private communication) 
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Rare earth ions exhibit CEF splittings which are usually very 
much smaller than similar splittings observed for ions of the iron 
transitions elements. In the iron transition series, the partially 
filled 3d electron shell is fully exposed to the CEF produced by 
surrounding ions, resulting in large CEF level splittings. The rela- 
tively small CEF level splittings observed for rare earth ions, 
which typically are of the order of a few hundred cm V probably 
arise because of large shielding effects resulting from the 5s°p^ 
electronic shells which surround the partially filled 4f shell. 

Present theoretical predictions of the influence of electronic 
shielding upon the CEF level splitting of rare earth electronic levels 
diverge. Burns (9) concluded that electronic shielding in the rare 
earth ions is of little importance rind that the difference between the 
CEF level splittings in the iron series and those in the rare earth 
series cannot be attributed to electronic shielding of the 4f electrons 
from the CEF by outer closed electron shells. In contrast, Eenander 
and Wong (10), Ray (11) and Watson and Freeman (12) conclude that 
electronic shielding plays a significant role in rare earth CEF level 
splittings . 

Quantitative estimates of actual shielding effects are hampered 
by the lack of sufficiently accurate atomic wave functions for rare 



(9) G. Burns, Phys. Rev. L2_8, 2121 ( 1962) 

(10) C.JL Eenander and E. Y. Wong, J . Chem . Phys . T8, 2750 ( 1963) 

(11) D.K. Ray, Proc. Phys. Soc. F2, 47 ( 1963) 

(12) R . E . Watson and A. J, Freeman, Phys . Rev. 13 3, A 1 57 1 ( 1964) 
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earth ions. Inadequate knowledge of the contributions of the core 
electrons is a primary source of uncertainty in our present under- 
standing of hyperfine interactions in rare earth (as well as in other) 
elements. Direct measurements of the influence of electronic 
shielding upon the nuclear hyperfine interactions and upon the CEF 
splittings of electronic levels therefore are highly desirable. 

This paper demonstrates the use of the technique of recoilless 
nuclear resonance absorption of gamma radiation as a means to ob- 
tain information on electronic shielding effects in rare earth isotopes. 
The procedure introduced here consists of combining measurements 
of the temperature dependent nuclear quadrupole interaction (per- 
formed by using the technique of recoilless resonance absorption) 
with measurements of the CEF level splittings (performed by using 
optical techniques). Specifically we report on determinations of the 
relevant electronic shielding factors for trivalcnt thulium based 

upon our gamma -absorption measurements of the nuclear quadrupole 

1 69 

interaction of Tm in thulium ethyl sulfate (13) and thulium oxide 
and on optical measurements of CEF levels by Wong and Richman 

(14), Gruber and Krupke (15), and Gruber et al. (16). 

1 £>9 

Tm appeared to be an isotope particularly suited for studies 

(13) A preliminary report of part of this work appeared elsewhere 
(7). 

(14) E . Y . W ong and I. Richman, J, Chem. Phys . 34 , 1 1 8Z (1961) 

(15) J.B. Gruber and W.F. Krupke, to be published 

(16) J. B. Gruber, W. F. Krupke and J. M. Poindexter, to be 
published 
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of electronic shielding, for the following major reasons: 

(1) The low energy of the 8.4 keV transition used results in a high 
Debye-Waller factor (recoil-free fraction) even at very high tempera- 
tures, thus permitting a measurement of the quadrupole interaction 
within an unusually wide temperature range. 

(Z) The separation of the excited levels belonging to the ground 

multiplet of thulium (L = 6; S - l) is rather large, with the first 

3 - 1 3 

excited level ( H^) some 5600 cm above the ground term ( H^). 

Thus the existence of the higher levels of the ground multiplet is of 
minor concern for the interpretation of our data in thulium, in con- 
trast to the situation prevailing in the case of some other rare earth 
ions . 

(3) The spin of the nuclear ground state (I = 1 /Z) and of the 8.4 keV 
excited state (I = 3/Z) is rather low , resulting in a small number of 
quadrupole hyperfine components of the gamma lines which are easily 
r es olvable . 

1 69 

(4) The nuclear collective model applies well to Tm thus per- 
mitting a rather reliable semi -theoretical estimate of the nuclear 

quadrupole moment of the 8.4 keV state. 

1 69 

(5) The relative abundance of Tm 



is 100%. 
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II. CRYSTAL ELECTRIC FIELD (CEF) INTERACTIONS 
A rare earth ion interacts in a salt with the CEF produced by all 
the ions which surround its position in the lattice. The dominant 
effect is the interaction of the CEF with the electrons in the partially 
filled 4f-shell. This interaction is weak compared to the spin-orbit 
interaction, in contrast with the situation prevailing in the case of 
iron-transition elements. As a result, the total angular momentum 
J remains a good quantum number for rare earth ions bound in 
crystals. The effect of the CEF then essentially is a partial or com- 
plete removal of the 2J + 1 fold spatial degeneracy of the orientation 
of which exists in a free ion. The actual number of electronic CEF 
levels depends on the symmetry of the field, while the level spacing 
depends on the strength of the interactions between the CEF and the 
4f electrons. The situation is illustrated in Fig. la. 

The potential energy describing the interaction between the CEF 
and a negative charge at position (r, $ , <p ) within the ion centered at 
the origin may be represented in good approximation by the following 
expansion, not including shielding from closed shells: 

+ n 

-eV(r,^,/)= y Y] A m r n tfP ('tf'.'f) , (D 

a — * * n n 

n m = — n 

if one assumes that there is no overlap between the charge distri- 
butions of different ions. In Eq. (1) the A™ represent lattice sums 
over point charges and effective multipole moments in the surround - 
ing ions. The relevant functions 0 m , which are linear combinations 
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Fig. la: Schematic of the atomic level splitting of a rare earth 

ion in the CEF. For a nuclear spin I = 3/2 the nuclear 
quadrupole interaction splits each CEF level into a 
doublet, which is the case illustrated. Typical overall 
CEF splittings are of the order of 10 ^ eV, while typical 
quadrupole hyperfme splittings are of the order of 10 ^ 
eV. 

Fig. lb: Schematic of the nuclear quadrupole splitting of the 

169 

8. 4 keV transition in Tm . The temperature dependent 
level splitting , which is typically of the order of 

10 ^ eV, is the average of the hyperfine splittings of 
Fig. la, weighted according to their Boltzman factors. 
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of spherical harmonics Y m and Y m , are defined as follows (17); 
r n n 



0^ =(2-4*6.. 

2n 



. . 2n)P, (cos'll 
2n 



* .m , , . (n— m) T ^m . n. 

<b = (2) mi 7 —r — Vr p cos • 

M (n + m) In 



(cos m^i 



n 



m 



> 0 



sin m 






where P and P m are Legendre polynomials and associated Legendre 



n n 



functions, respectively. In particular, we obtain for n = 2: 



0 ^ - 3 cos z ^ - 1 



(2a) 



0^ = sin a if cos 2 



(Zb) 



0 sin 1 ^ sin 2 



(2c) 



Specifically, the Hamiltonian describing the interaction between 
the CEF and the electrons in the partially filled 4f shell of rare earth 
ions, including the effect of shielding via the closed electron shells 
of the central ion is given by: 



H 



(4 f) 
CEF 



E 

k 



£ A" 

n m 




S 

n 







( 3 ) 



(17) The normalization of the functions 0 m ( 7^, V 7 ) is arbitrary; 
the choice adopted here is the one most commonly used 



The terms proportional to r™ describe the potential energy due to 
the direct interaction of the CEP" with the k-th electron in the 4f-shell 
while the terms proportional to S^(r^) describe the additional poten- 
tial energy arising from a deformation of the closed electron shells. 

The interaction described by the Hamiltonian in Eq.(3) splits the 
electronic ground state of the free ion, characterized by total angu- 
lar momentum J , into a number of CEF levels. We shall assume 
in calculating these CEF levels that the angular and radial parts of 
the free ion wave functions can be factorized and that higher terms 
with diffe rent values can be neglected. Under these circumstances 
we are dealing with a manifold of states belonging to the same J and 
it is then convenient to replace the angular operators occurring in 
the Hamiltonian, Eq.(3), by equivalent operators (18). The relevant 
matrix elements then are of the form 



H = V A m /r n >„ <J || d || J> <J,m T |o m (J , J , J (4) 

m T m Tt Ls n \ /E x 11 n" ' \ J‘ n x y z 1 J 
J J' n m 

where ^/^ = ( 1 -%) 4"^ < 5 ) 



(18) K.W.H. Stevens, Proc. Phys. Soc. Ao5, 209(195Z); 

R.J. Elliott and K. W. H. Stevens, Proc. Roy. Soc . AZ 1 8 , 
553 (1953); J.P. Elliott, 13. R. Judd and W. A. Runciman , 
ibid. AZ40 , 509 ( 1957); R. Orbach, ibid. AZ64 , 458(1961) 
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and 



O' 

n 


< U 4fl 


s „< r > kr >/<' r " 


>4 f 


(6) 




< U 4f 


1 r 1 U.k 

1 nl 4f/ 




(7) 


II 


a » ft, 


y for n = 2 , 4 , 


6 





In these expressions is the radial part of the electronic wave 



functions for the 4f-shell. The functions O m (J , J , J ) are opera- 

n x y z ^ 

tor equivalents; those relevant for this work are listed in Table I. 
The expressions <^T |j / Z^ ||j are reduced matrix elements (19), 
which for the more general case of intermediate coupling are avail- 
able for Tm^ in the literature ( 14)-( 17), (20). 

It is in principle possible to calculate the parameters A™ and 
but difficult in practice. Difficulties are in the evaluation 
of the ’'lattice sums" A™ because of a lack of sufficient knowledge of 
the ionic position coordinates and their temperature dependence as 
well as of the values of moments in the surrounding ions (21). The 
evaluation of the radial integrals e> E . which are the expectation 
values of r n for the 4f shell modified by contributions from closed 
shells to the electric multipole fields at the 4f electron positions, is 
hampered by the lack of knowledge of sufficiently accurate atomic 



(19) B.R. Judd, Proc. Roy. Soc. A24 1 , 414 ( 1957) 

(20) J. B. Gruber and J.G. Conway, J. Chem. Phys. 32 , 1 53 1 
( I960) 

(21) M.T. Hutchings and D. K. Ray, Proc. Phys. Soc. 81,663 
(1963) 
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- 13 - 



fM N 

‘"’I 



m 

I 



o 

fM 



S3 

r- 



4 > 

I 

<m n 
4 



m 

fM 

lT> 



4 

►— > 



m 

4 



'“l 

xO 



M 



LD 

fM 



fM I 



y 

4- 



fM I 



4>\ 



l r 

1 



fM I 

I 






fM I 

I 

<vj 4 



m 

1 



fM 

4 > 

Ln 

o 



tT N 4 

» — 1 1 O 

vO 

fM 

I CT 







fM I 


4 


4 


fM 1 


1 


o 

4 






fM 4- 


, 


4 

4) 


fM + 




1 








4 

z 


1 





• H 

fM 




fM SI 


fM 




lO 








rn 


4 




i— < 




II 


II 


II 


11 


O fM 


fM fM 


(M 

1 fM 


O Tf 


0 


o 


0 


o 



I 

fM S3 



(M 

o 



4 

I 

fM NJ 
4 



fM 

i 



4 

4 4 



o 



4 4< 

o 



I 

rf 4 

m 

fM 



I rT 

0 



in 

o 



vO S3 

> 4 * 



m 

fM 



O vD 

o 



LT) 

I 



(continued) 



TABLE I. Operator Equivalents (continued) 



14 



I 



rvj + 



4 ] 

O 



ro 

O 



«M I 

I 



4 



4 ] 

O 



rsj 

O 



oo 



oo 

rO 



*-5 

O 



>“5 



oo 



I 

4 S3 

rO 



4 



4 
*— ) 



> 

o 



4 

4> 



I 

4 S3 

ro 

rO 



eg l 
•“>* 

4 

<vj 4 - 



4 
*— ) 



o 

4 



I 

4> 



*— > 

00 



I 

4 N 

m 

rO 



4 

I 

<vj S3 






I 

4 N 

rO 

rO 



rg I 

■"M 

I 

4 



4 

I 

rg n 



4 



4 

I 

rvj s; 



I 

4 4 



4 


4 


eg 


1 j 

4 


4 4 
4\ 


4 


4 


^ S3 


4 


eg N 






ro tSJ 


ro 

4) 


4 

<VJ S3 


co 

4 


4 

4~ i 


4 ' 

1 

4 4 


m 


4 


rO 


4 


4- 


f ) 


r\J 


- — - 


41 





4 

+ 






’ — * 




. — i 


2d 


* 


4 


4 


4 


4 


r > 




4 




4) 




00 




s ‘ w " 






, — , 


co 


-f 


4 


4 


) 


00 


1 




00 


4 


00 


rO 


_ _ 


r J 


1 4 




— 


1 





eg S3 



vO I 



vO ‘ 


1 


4 


vO 4 


vO 4 




^1 


, — , 




« r4 

41 


41 






«— H 



M 



II 



II II II 



II 



4j vO 

o 



4J 

• nO 



o 



4 vO 

o 



4 

I vO 



0 



vO o 

O 



vO 
I 4D 



O 



-15- 



wave functions for bound rare earth ions. For these reasons it is 
therefore preferable to introduce the "CEF parameters" 



C m = A 
n 



m / n\ 
n V / J 



( 8 ) 



to be determined by experiment. The point symmetry of the central 
ion drastically reduces the number of CEF parameters (ZZ). In the 
case of rare earth ions only the terms with n = Z, 4, 6 need to be 
considered, with the effects of n = 1 , 3, 5 being negligible in most 
cases ( ZZ) . 

The wave functions of the V'-th CEF level will be taken as a 

linear combination of eigenvectors of the total angular momentum J , 



Vy* U 4r (r) E a - 

(m j) 



(m T ) , m 
J 



J IjJ 



(9) 



The expansion coefficients and the energy eigenvalues E ^ 

follow from the diagonalization of the interaction matrix H 
Eq. (4). 






(ZZ) A compilation of the relevant values n and m for various 
crystal symmetries was given by J. E. Prather, NBS 
Monograph 19 (1961) 
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III. THE NUCLEAR QUADRUPOLE INTERACTION 
Each of the CEF levels may produce a magnetic field and an 
electric field gradient at the nuclear site; this results in hyperfine 
splittings of the electronic levels. A rare earth nucleus thus ex- 
periences at a certain time a magnetic field and an electric field 
gradient which depends on the electronic state that is actually popu- 
lated at this time. The situation substantially simplifies at elevated 
temperatures where the spin-lattice relaxation phenomenon produces 
rapid transition between the different CEF levels. The nucleus 
under these circumstances experiences a magnetic field and an elec- 
tric field gradient which in essence result from averaging these 
fields over all electronic states weighted according to the population 
numbers. This averaging process, which essentially constitutes a 
time averaging process, holds only if the significant electron relax- 
ation times are short compared to all other relevant times sucli as 
the nuclear lifetimes and the nuclear precession times, a situation 
prevailing at temperatures above a few degrees Kelvin. In particu- 
lar, the magnetic hyperfine interaction cancels in the absence of an 
external magnetic field and all one is left with is the quadrupole 
hyperfine interaction (23), (24). An example of this situation is 
illustrated in Fig. lb for an assembly of nuclei. The quadrupole 

(23) R. L. Cohen, U. Hauser and R. L. Mdssbauer, Pr oc . 
MHssbauer Conf. 2nd , (John Wiley and Sons, N. Y. , 1962) 
p. 1 72 

(24) R. L. MOssbauer, Rev. Mod. Phys. 36^ 362(1964) 
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interaction is strongly temperature dependent since the overall CEF 
splitting within the lowest electronic term is only of the order of a 
few hundred cm ^ . 

The electric field gradient which interacts with the nuclear 
quadrupole moment of a rare earth nucleus bound in an ionic crystal 
has four significant sources: 

1) One contribution is the direct field gradient produced at the 

nuclear site by all of the ions surrounding the host ion which 
contains the nucleus in question. This contribution in the 
case of rare earth ions is usually negligible in comparison 
with the contributions from other sources, particularly at 
low temperatures. 

Z) Another contribution results from the electric field gradient 

produced at the nuclear site by the electrons in the partially 
filled 4f-shell of the host ion. This field gradient results 
from the interaction of the 4f-electrons with the CEF pro- 
duced by the surrounding ions. This interaction eflectively 
induces electric mu lti pole moments (multipole polarization) 
in the 4f-shell; the quadrupole part of this polarization con- 
tributes to the electric lield gradient experienced by the 
nucleus. 

3) A distortion is usually also induced by the CEF in the closed 

electron shells, yielding another contribution to the total 
field gradient experienced by the nucleus. This contribution 
is proportional to source (1), with proportionality factor -y^. 
The absolute value of the proportionality factor is in the 
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case of the rare earths usually large in comparison with 
unity, thereby leading to such an enhanced field gradient 
(antishielding effect) that it often becomes comparable with 
the one resulting from source 2). This is the "lattice" 
Sternheimer effect (2S)-(28). 

4) Another field gradient contribution due to an induced quadru- 

pole moment in the closed electron shells results from the 
interaction of the closed electron shells with the electrons in 
the partially filled 4f-shell. This relatively small contri- 
bution, which is proportional to source 2), with proportionali- 
ty factor -Rq, is the "atomic" Sternheimer effect (5), (28). 

Collecting the different contributions, we obtain for any compo- 
nent ea . of the electric field gradient tensor 



eq . 



= u 



-y ) eq (Lat) 

GO 1J 



<‘-V ^ 



i, j = 1 , 2, 3 



( 10 ) 



where y and R„ are the lattice and atomic Sternheimer factors, 
' co Q 

respectively, as introduced above. 

(25) E.G. Wikner and G. Burns, Phys. Letters 2 , 225(1962) 

(26) D.K. Ray, Proc. Phys. Soc._82, 47 (1963) 

(27) R.M. Sternheimer, Phys. Rev. 132 , 1637 (1963) 



(28) A.J. Freeman and R.E. Watson, Phys. Rev. 132 , 706 
(1963) 
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In the principle axes system of the electric field gradient tensor 

(V ) 

the nuclear quadrupole interaction Hamiltonian Hq associated with 
the (V)-th CEF level of the ion is given by 

(ID 

(V)_ e z Q 



H 



Q 



41(21 



17 , ( [< 1 - r q> 






where 1^, JL , are the usual nuclear spin operators and Q is the 

nuclear quadrupole moment. The quantities and 

determine the direct contributions to the electric field gradient at 
the nuclear site produced by the surrounding ions in the lattice and 
by the 4f-electrons of the host ion, respectively. They are defined 
by 

eq { .^ at) = [c> 2 V(r,^, </> ) / 3x. d x. ] f = Q , (12) 



where V(r , ^ , V* ) is defined in Eq. (1), and 

4f -elect rons 



13) 



= < y l5ZZZIZ[5 2 ( ij - "r - | ) / 3x i 9 x iI r --o| y ) 



where the wave function | of the \/-th CEF level is of the form 
given by Eq. (9) . 

Explicitly we obtain for the lattice contribution from Eqs.(l), 

( 2 ), ( 12 ) 



2 (Lat) ..0 z, (Lat) (Eat), . A 2 

e 2 q' -4A ; e 2 (q' - q' ') = - 4A 

n zz 2 xx yy ' 2 



(14) 
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Expressing the contributions from the 4f -electrons within a manifold 
of states of constant ^ in terms of operator equivalents, we obtain 
from Eq. (13) 

<Hs ( z 4 z f) |^ =-< J Ik II J > < r ' 3 >4f <^l 3 4 2 -JL 2 K) 

<H-9x 4 x f) -5yy f V> = - (3/2) < j ll-ll J > < r " 3 >4f <^lif 

(15b) 

where <^r is defined by Eq.(7). 

Usually one observes only an average field gradient from the 4f- 
electrons, which is a field gradient from the individual CEF levels 
weighted according to their Boltzman factors, as discussed above. 

If we consider only those electronic states which belong to the lowest 
manifold spanned by the state vector , then the average direct 
contribution from the 4f-electrons to the electric field gradient acting 
on the nucleus at temperature T is given by (29) 

2 J + 1 

JC <Hs! 4f V> ‘ ex P ( “E i/ /kT) 

- — ( 16 ) 

2 J + 1 

FI exp(-E y / kT) 

V = 1 

The diagonal component of the averaged total electric field gradient 
tensor is according to Eqs.(10) and ( 1 6 ) given by 



(4f) 

<A« > 



(29) A more general description including effects of higher J 
states is given in the Appendix I 
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/i , ( Lat) 

( 1 - y ) q:. 

oo 11 



“-V <s‘ i ,t) >T 



(17) 



where we have neglected any temperature dependence of the lattice 

* * K ( Lat ) 

contribution q. . 

H 

The total Hamiltonian describing the average quadrupole inter- 
action at temperature T may now according to Eqs. ( 1 1) and (17) be 
written as 

( 18 ) 



H q (T) = 



c z Q 



41(21-1) 



<a zz > T < 3 i- 2 z -l 2 > + <a X x-a yy >T Hi*** 



We shall now apply the preceding formalism to the particular 

1 69 

case of Tm . The twofold degeneracy of the nuclear ground state 

1 69 

of Tm (I = 1 /Z) is not removed by the Hamiltonian, Eq. (18); the 
8.4 keV excited state (I = 3/Z), on the other hand, is split by the 
nuclear quadrupole interaction into two states. Their energy separa- 
tion ^AE^> ^ which follows from the diagonalization of the Hamiltonian 
Hq(T) is given by 




-f 



3 




l 

z 



(19) 



This expression may be written in more detail, by using Eqs. (5), (8), 
(14), (15), (17): 



-LZ- 



(AE^ = (e 2 Q/2){ [<J||a|J> <r % 



<3J 2 - J 2 > 
z I 



4C 



<^ 2 >E 



( 1 - y )] 

ou i 



M-.j] )' 

UO) 

where 3 J 2 - J 2 /L and <^J 2 + J 2 Xr- are thermal averages defined 

~-z ~~ T — T 

as those given by Eq.(l6), while the parameter (r ^ is defined by 



a 



3 

Z 



<J|| a\\ J> < r " 3 > Q <J 2 + + J 2 _ > T + 



4c: 



2 <r 2 \- 



<r' 3 > Q = (1- R q ) <r " 3 > 



4f 
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It is just this splitting (AE) r ^ that is measured as a separation of 
gamma lines in recoilless resonance absorption experiments. 

Several additional hyperfine interaction mechanisms which 
contribute to the net nuclear quadrupole coupling of a rare earth ion 
have been neglected in our calculations. These additional contri- 
butions arise in second-order perturbation theory with the principal 
effects coming from the magnetic hyperfine interaction itself (30) 
(the so-called pseudo -quadrupole coupling) and from the quadrupole 
interaction with states of higher admixed into the ground state 
multiplet by the CEF. We have made calculations of these contri- 
butions for the compounds covered in this paper and they amount to 
less than 1 % of the total quadrupole interaction energy. 

In order to compare experimental results with theory within the 
framework of the CEF model it is convenient to replace in the theo- 



(30) R.J. Elliott, Proc. Phys. Soc. B70 , 1 19 ( 1957) 
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retical expression for the quadrupoie splitting ( AE)^ all quantities 
involving electronic radial integrals (the theoretical determinations 
of which is presently somewhat uncertain) as well as the nuclear 
quadrupoie moment by experimentally observable parameters* For 
this purpose we introduce the dimensionless parameters 



P, = e*Q <r' 3 > Q <J |a |j > /C° 


(22a) 




(22b) 



Expressed in terms of these parameters the quadrupoie splitting 

1 69 

in Tm reduces to 



<AE> t 




p l 



< 3 J 2 -J 2 > x + 4C 



2 



+ 



(23) 



+ 



3 



(3/2)C° p 



< J + 



+ J 2 >, 



- T 



+ 4C 




2 



The temperature averages ^ anc ^ + 

the framework of the CEF model depend only on the experimentally 

observable CEF parameters C m , 

* n 
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IV EXPERIMENTAL TECHNIQUE 



The nuclear quadrupole interaction was measured by using the 

technique of the recoilless nuclear resonance absorption of gamma 

1 69 

radiation (2). The partial decay scheme of Er is shown in Fig. 2. 
Measurements of the gamma resonance absorption were performed 
as a function of the relative velocity between sources and absorbers. 
The measurements involved sources of erbium trifluoride (ErF^) and 
erbium oxide (Er^O^) and absorbers of thulium ethyl sulfate 
(Tm (C^H^SO^)^ 9H^O . abbreviated to TmES)and thulium oxide 



Anhydrous ErF^ provides an excellent source for experiments 

1 69 

utilizing the 8.4 keV line of Tm . The crystal structure of the 
heavy rare-earth trifluorides has been investigated by Zalkin and 
Templeton (3 1). At temperatures below about 900 — 1000° C the 
stable phase is orthorhombic, space group -Pnma, having four 

formula units per unit cell. The rare earth ions are c r ystallog ra - 
phically equivalent, having the point symmetry m. Thus, although 
the electric field gradient tensor eq . . is not axially symmetric, all 
erbium (or thulium) nuclei experience the same ec^... Therefore, 
the quadrupole splitting of the recoilless absorption line given by 
Eq.(23), may be expected to pass through zero or at least through a 
minimum at a specific temperature (550° K in this case). The advan- 
tages of a single-line source are thereby obtained. The line width 
obtained this way with sources of Erf' ^ is less than with sources of 



(31) A. Zalkin and D. H. Templeton, J. Am. Chem. Soc. 7 5 , 
2453 (1953) 
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Fig. 2: 



Partial decay scheme of Er 
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where the presence of non -equi valent erbium sites complicates 
the situation (23), (24). At the same time, with reasonable precau- 
tions, the ErF^ can be maintained at the critical temperature for 
periods of several weeks without decomposition or reaction. This 
chemical stability does not exist with most other erbium salts in 
which the erbium ions are also c rystallographically equivalent (e. g. , 
the sulfate, nitrate, chloride). 

Anhydrous ErF-j was prepared from erbium metal or erbium 
oxide by a "wet” process. The metal or oxide was first dissolved in 
a small quantity of nitric or hydrochloric acid in polyethylene centri- 
fuge tube. A few ml of aqueous hydrofluoric acid were then added 
and the mixture heated at approximately 100° C in a water bath for 30 
minutes. The somewhat gelatinous ErF^ precipitate was then centri- 
fuged down, the excess solution decanted off, the precipitate washed 
with distilled water, centrifuged three to five times and dried in air 
at roughly 100° C. Air drying yields a hydrated ErF^ of unknown 
composition. To remove the water of hydration, the dry contents of 
the centrifuge tube bottom were transferred to a small tantalum boat 
and annealed in an evacuated fused quartz tube. Experience showed 
that the hydrated ErF^ could be converted directly into a mixture of 
the several forms of oxyfluoride (32) if the annealing temperature 
was raised too rapidly. The procedure finally adopted was to hold 
the hydrate at room temperature at about 10 5 torr for at least 12 
hours in order to pump off most of the water. The temperature was 



(32) W. H. Zachariasen, Acta Cryst. 4, 231 ( 1951) 
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then raised slowly (in 6 hours) to 150°C, thus removing virtually all 
of the water. Finally, the temperature was raised to 850° C in an- 
other 6 hours and then reduced back to room temperature within 2 
hours. This procedure yielded consistently good clean x-ray powder 
patterns of the orthorhombic phase without a trace of the hexagonal 
phase appearing (31). ErF^ prepared in this manner appears to re- 
main stable at room temperature over an indefinite period of time. 

At elevated temperatures care must be ex/^ercised to avoid reaction 

with oxygen or water vapor. Sources of ErF^ were prepared in the 

1 1 > 8 

above manner from Er^O^ (usually enriched in Er ) or from erbium 
metal after irradiation in the Materials Testing Reactor, Arco, Idaho. 
Alternatively, the ErF^ was prepared first and then irradiated. Iden- 
tical spectra were obtained by the two methods. 

Absorbers of TmF^ were used in order to experimentally determine 
the critical temperature at which the narrowest possible emission 
line is obtained with sources of ErF^. Figure 3 shows the tempera- 
ture dependence of the quadrupole splitting in TmF^. The source was 
mounted in a small evacuated oven shown in Fig. 4. The absorber 
was maintained in a helium atmosphere within an oven equipped with 
beryllium windows. 

It is interesting to note that the same minimum line width 
(1.8 cm / sec) was obtained in both the trifluoride - trifluoride and tri - 
fluoride - ethyls ulf ate experiments. This strongly suggests that the 
quadrupole splitting of the trifluoride line does indeed pass very near 
to zero at 550° K (24). This minimum observed line width of 
1. 8 cm/sec may be compared with the theoretically predicted line 



-28- 



Fig. 3: 



Temperature dependence of the quadrupolc splitting 

1 69 

of the 8.4 keV gamma line of Tm using an ErF^ 
source and a TmF^ absorber. Source and ^lbsorber 
were maintained at the same temperature. 



QUADRUPOLE SPLITTING <AE> 
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Fig. 4: 



Details of moveable source oven. The entire oven 
(weight 0. Zkg) was moved relative to the absorber 
by the cam drive. The main body of the oven was 
made of stainless steel. A similar oven made of 
aluminum was used with the transducer drive. Tin* 
heating element was fabricated from nichrome strips, 
1/16 x 0.005 in. For a source temperature of 5 50° K 
the power input was 30 W. 
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width. In a t r ans mis s ion - vs - veloc ity measurement since an emission 
line of width T is moved over an absorption line of width V , one ex- 
pects a minimum line width of 2T . In the case at hand this is (in 
velocity units) 



2 T c/E - Z'fic/'tE, - 0,74 cm/sec 

y y 

based on the lifetime ( 3 3) , ^ = 6.28* 10 ^ sec, and energy, 

169 

E “ 8.42 keV, of the first excited state of Tm . However for an 

y 

absorber of finite thickness this width increases by a correction 
factor which is 1,47 in our case for an absorber of 5 mg/cm z of 
thulium and a total conversion coefficient of 325 from Kankeleit et al. 
(34). Thus the predicted line width is 1.09 cm/sec which must be 
compared with the minimum observed line width of 1.8 cm/sec. The 
observed line width is 1. 6 times broader than expected. The origin 
of this line broadening is uncertain. 

Absorbers of TmES were prepared by crushing single crystals. 
Absorbers of Tm^O^ and sources of (enriched) Er^O^ were prepared 
from commercially available material. Absorbers of all materials 
to be used below room temperature were prepared by mixing the 
powdered samples with a soft wax and pressing the mixture into thin 
disks between mylar films. Absorbers and sources of all materials 

(33) R. E. McAdams, G. W. Eakins, E.N. Hatch, Phys. Letters 
6, 219 ( 1963) 

(34) E. Kankeleit, F. Boehm, R. Hager, Phys. Rev. in press 
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to be used above room temperature were prepared by settling the 
powdered samples from a slurry of dry acetone onto l/Z mm thick 
beryllium windows. 

The relative velocities required for Dopple r - sh i fti ng the gamma 
lines were produced by using both cam drives (35) (providing constant 
velocities) and transducer drives ( 36) (providi ng constant acceleration). 
A block diagram showing the experimental apparatus for use with the 
cam drive is shown in Fig. 5. The experimental arrangement used 
with the transducer drive is given in Ref. (36). Proportional counters 
filled with one atmosphere of a mixture of 90% argon and 10% methane 
(by volume) and equipped with l/Z mm thick beryllium windows were 
used as detectors, see Fig. 7. 

A cryostat specifically designed for recoilless resonance absorp- 
tion experiments with low energy gamma radiation was used for the 
measurements (37). The sample temperatures in the range from 1 0°K 
up to 300° K were attained by either controlled heating of the cooled 
sample holder, by pumping on liquified gases, or by using exchange 
gas cooling. The sample disks were clamped between thin beryllium 
disks soldered to the cryostat sample holder in order to ensure good 
temperature uniformity and stability. Temperature measurements 
were made using carbon resistors and thermocouples. The oven used 

for heating sources to Z000° K is shown in I ig. 8. 

1 69 

Some typical Mbssbauer spectra for Tin are shown in Fig. 9. 

(35) R.L. Mbssbauer, Proc. Mbssbauer Conf. Znd , (John Wiley 
and Sons , New York^ 1 9 6 Z ) , pTift 

(36) E. Kankeleit , Rev. Sci. Instr. 35_, 194(1964) 

(37) F. T. Snively, to be published 
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Fig. 5: 



Block diagram of experimental apparatus for use with 
the cam drive. A detailed schematic of the pro- 
grammer is presented in Fig. 6. All other electronic 
equipment shown is commercially available. The cam 
drive has been described elsewhere (35). The recoil- 
less resonance absorption measurements were per- 
formed by first moving the source at a normalizing 
speed (12 cm /sec), then at a measuring speed (v), and 
finally at the normalizing speed. The three runs were 
of about 5 min each. The counting rates, C, during 
the two normalizing runs were averaged and combined 
with the results of the measuring speed to yield the 
amount of absorption, A(v). 

A(±v) = [ C ( ± 1Z) AVE -C(±v)j /C(± 12 ) ave 
T his sequence of events was repeated until the counting 
statistics were satisfactory. The information at the end 
of each run was printed out on a typewriter for moni- 
toring purposes and punched out on paper tape for 
processing by an electronic computer. The entire 
process was automatic except for changing from one 
measuring speed to another. 
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Fig. 6: 



Schematic of programmer for use with cam drive. 
This programmer controlled the changes from the 
normalizing speed to the measuring speeds and 
initiated the counting and print-out cycles of the 
scaler, see Fig. 5. 
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FIG. 6 
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Fig. 7: 



Proportional counter. The anode of the counter is 
a 3 mil stainless steel wire kept under light tension 
by a spring in the end fitting. The beryllium window 
is 0. 5 mm thick. This counter was operated at about 
2500 V with 1 atm of 90% argon and 10% methane, by 
volume. Under these conditions the resolution at 
8.42 keV was about 17%. 
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Fig. 8: 



High temperature, resistance heating, source oven. 

The heating coil was fabricated from two pieces of 

56 mil tungsten wire 15 in long. These two pieces 

were connected in parallel and the source container 

was suspended from the midpoint of each piece. For 

a source temperature of 1970° K the power input was 

-4 

4. 8 kW. A vacuum of 10 torr was maintained in 
the oven. The thermocouple made of tungsten vs. 
tungsten - 26% rhenium is reliable to 3100° K. 
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Fig. 9: 
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Quadrupole splitting of the 8.4 keV level of Tm 

in an absorber of thulium ethyl sulfate (5 mg/cm 2 

1 69 

of thulium). A ’’single line” source of Er in ErF^ 
was used at the critical temperature T = 550° K 
throughout curves a-d. The spectra a, b and c, d 
were obtained by using a cam drive and a transducer 
drive, respectively. 
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FIG. 9. 



